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We describe the design and current performance of the high-flux backscattering spectrometer 
located at the NIST Center for Neutron Research. The design incorporates several state-of-the-art 
neutron optical devices to achieve the highest flux on sample possible while maintaining an energy 
resolution of less than 1 fieV. Foremost among these is a novel phase-space transformation chopper 
that significantly reduces the mismatch between the beam divergences of the primary and secondary 
parts of the instrument. This resolves a long-standing problem of backscattering spectrometers, and 
produces a relative gain in neutron flux of 4.2. A high-speed Doppler-driven monochromator system 
has been built that is capable of achieving energy transfers of up to ±50 fieV, thereby extending 
the dynamic range of this type of spectrometer by more than a factor of two over that of other 
reactor-based backscattering instruments. 

PACS numbers: 61.12.Ex.61.12.-q 
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I. INTRODUCTION 



Neutron scattering is an invaluable tool for studies of 
the structural and dynamical properties of condensed 
matter. Neutron sources produce neutrons with wave- 
lengths that span the interatomic spacings in solids or 
the diameter of complex macromolecules, while at the 
same time having energies that match, respectively, the 
lattice vibrational frequencies in solids or the slow diffu- 
sive motions of atoms. The particular neutron scattering 
technique known as backscattering Jl], |) is able to resolve 
energies below 1 [leV, which is well beyond the reach of 
conventional triple-axis and neutron time-of-flight spec- 
trometers. Thus neutron backscattering spectroscopy is 
ideally suited to the study of dynamics such as slow mo- 
tions in complex liquids, jump diffusion, and quantum 
rotational tunneling. 

The principle limitation of backscattering spectrom- 
eters has long been the relatively low neutron flux on 
sample that they produce. This is, of course, a direct 
consequence of the excellent energy resolution they pro- 
vide. In this paper we report on the design and per- 
formance of the new high-flux backscattering spectrom- 
eter (HFBS) located at the NIST Center for Neutron 
Research that addresses this limitation. Compared to 
other backscattering spectrometers, the HFBS delivers a 
higher neutron flux to the sample in large part by the 
use of a novel de vice ca lled a phase-space transforma- 
tion chopper (Sec. IIIC). In addition, a newly designed 
Doppler-driven monochromator (Sec. [HE), which oper- 
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ates with a cam machined to produce a triangular veloc- 
ity profile, extends the dynamic range of the spectrom- 
eter by more than a factor of two beyond that of other 
similar instruments. In a departure from other backscat- 
tering spectrometers, the scattering chamber is operated 
under vacuum instead of in an argon or helium gas envi- 
ronment, which improves the signal-to-background ratio 
substantially. 



II. BACKSCATTERING 

Backscattering spectroscopy exploits the fact that the 
wavelength spread AA of a Bragg-diffracted neutron 
beam decreases as the scattering angle 28 approaches 
180° (see Fig. |l|). This is easily shown |J by differen- 
tiating Bragg's law (A = 2c? sin 6), and then dividing the 
result by A to obtain 
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As 8 — > 90° the angular term vanishes. This results in 
a minimum value of AA/A, and hence the energy resolu- 
tion, that depends on the spread Ad and average value 
d of the lattice spacing between crystal Bragg planes. In 
the kinematic limit this minimum is zero. However, dy- 
namical scattering theory shows that the lattice gradient 
term Ad/d is non-zero, even for perfect single crystals. In 
this case, Ad/d is given by the Darwin width of the reflec- 
tion being used to monochromate the neutron beam. Q] 
This presents a fundamental lower bound for the energy 
resolution that can be obtained via backscattering, which 
depends entirely on the structure factor of the reflection 
being used to monochromate the beam, and the number 
density of unit cells within the monochromating material. 
In the backscattering condition the neutron beam is nor- 
mal to the Bragg planes, corresponding to a Bragg angle 
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FIG. 1: Illustration of the backscattering principle. Verti- 
cal shaded regions correspond to equal angular spreads A8, 
but vastly different wavelength spreads AA depending on the 
Bragg angle 9. As 6 — > 90°, AA (and thus the energy resolu- 
tion) approaches a minimum. 

9 = 90°. However, the neutron trajectories in a beam 
are never perfectly parallel. Therefore, some neutrons 
will strike the crystal Bragg planes at angles slightly 
less than 90°, thereby satisfying the Bragg condition at 
different values of A. Consequently the spread A9 in 
incident angle will also contribute to AA. p{ Note that A9 
isn't necessarily (and usually is not) equal to the beam 
divergence, as it is most often set by the ratio of the 
source size to the distance between source and Bragg 
planes. If A9, is small, then 

AE AA /Ad 1, A ,o\ 

Most backscattering instruments use the {111} lattice 
planes of perfect silicon crystals to monochromate the 
incident beam as well as to analyze the energy of the 
scattered beam. This is true for the HFBS as well, 
so for the sake of convenience we define Xq = 2d — 
6.2712A, @ k = 2tt/A = l.OOA" 1 , v = 630.8m/s, 
and Eq — 2.08 meV. In this case the lattice gradient term 
Ad/d = 1.86 x 10~ 5 . As an example of how much the 
angular spread A9 contributes to the energy resolution, 
one would need a A9 = 0.70° to match the lattice gra- 
dient contribution to the energy resolution, which is a 
small angular spread for a neutron beam. Equation (|^) 
would then imply an energy resolution of about 0.16 fieV 
for the diffracted beam. 

The maximum momentum transfer accessible given 
neutrons of wavelength Ao is Q = 47r/Ao = 2.00A _1 , 
whereas practical considerations generally limit the min- 
imum useful Q to ~ 0.1 A^ 1 due to the non-zero diver- 
gence of the neutron beam incident on the sample. The 
energy range over which the dynamical properties of a 
sample can be studied is set by how much the energies of 
the incident and scattered neutron beams can be shifted 
relative to each other. This shift cannot be achieved 



by varying the Bragg angle of the monochromator, as 
is often done on a triple-axis spectrometer, because do- 
ing so ruins the excellent energy resolution. Instead, en- 
ergy transfers are obtained using other methods such as 
varying the temperature of the monochromator with re- 
spect to that of the analyzer, resulting in a continuous 
change in the d-spacing, or via a Doppler motion of 
the monochromator crystal, Q which is the method cho- 
sen for the HFBS. Since the analyzer crystals are fixed, a 
backscattering spectrometer can be compared to a triple- 
axis instrument operating in a fixed final energy config- 
uration. This is also often referred to as an inverted 
geometry configuration. Typical backscattering instru- 
ments with sub-yLteV resolution can reach energy transfers 
from ±10 to ±15 /xeV using earlier style Doppler-driven 
monochromator systems. [0 



III. GENERAL SPECTROMETER LAYOUT 

The design of the HFBS backscattering spectrometer is 
optimized to provide a large dynamic range and the high- 
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FIG. 2: General layout of the NIST Center for Neutron 
Research high-flux backscattering spectrometer (HFBS). 
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est neutron flux on sample possible while maintaining a 
sub-/KeV energy resolution. |3] To achieve these goals the 
HFBS design incorporates several state-of-the-art neu- 
tron optic devices, which are identified in the schematic 
diagram of the spectrometer shown in Fig. ^. Neutrons 
from the cold source of the 20 MW NCNR research re- 



actor (Sec. [II A) are conducted along a 41.1m straight 
neutron guide that is 15 cm high by 6 cm wide, and pass 
through beryllium and bismuth filters and a velocity se- 
lector. A converging guide (Sec. IIIB ), located after the 
local beam shutter, focuses the neutron beam cross sec- 
tion down to 2.8 cm x 2.8 cm, which enhances the neu- 
tron flux by ~ 3.9. The neutrons then enc ounte r a phase 
space transformation (PST) chopper (Sec. Ill C ), a device 
that Doppler-shifts the incident neutron wavelength dis- 
tribution towards the desired backscattered wavelength 
Ao- The PST chopper provides an additional gain of 4.2 
in neutron flux, but at the expense of a sizable increase 
in divergence. However in doing so the PST chopper 
alleviates the severe divergence mismatch between the 
primary and secondary spectrometers, which has been a 
long-standing problem with backscattering instruments. 

Continuing on from the PST chopper, the neutrons are 
backscattered from a spherically focusing monochroma- 
tor, strike the sample, and then are backscattered a sec- 
ond time from a spherically focusing analyzer system be- 
fore they finally reach the detectors. Note that backscat- 
tered neutrons must pass throug h the sample twice. Both 
the monochromator (Sec. HID) and analyzer (Sec. [II F) 
are composed of large, bent, silicon {111} crystals. These 
crystals are intentionally bent to increase the lattice gra- 
dient term in the energy resolution (see Eq. (|^)) in order 
to obtain a roughly three-fold increase in neutron count 
rate, while maintaining the sub-/ieV energy resolution 
constraint. Finally, a large energy difference (dynamic 
range) of up to ±50 /^eV can be established between the 
monochromated and analyzed beams using a cam-based 
Doppler drive that pr oduce s an oscillatory motion of the 
monochromator (Sec. HIE). 

The excellent energy resolution of backscattering in- 
struments comes at the cost of an inherently low neu- 
tron flux on the sample. But gains provided by the con- 
verging guide, the PST chopper, and the strained silicon 
monochromator and analyzer crystals help to compen- 
sate for this. In addition, the HFBS analyzer stands 
2 m tall, spans 165° in 29, and subtends nearly 23% of 
4-7T steradians (Sec. Ill F ) , making it the largest analyzer 
of any other backscattering instrument. It is composed 
of about 12 m 2 of silicon. This large analyzer provides 
yet another gain over other spectrometers through an 
increased count rate that is not included by measure- 
ments of the neutron flux on the sample. In the follow- 
ing subsections wc discuss the various components of the 
spectrometer in detail starting from the cold source and 
ending with the detectors. In Sec. IV we report on the 
performance of the HFBS and present results from sev- 
eral experiments. Recent upgrades and some potential 
prospects for improvements to the spectrometer are dis- 



cussed last in Sec. [v| 



A. Neutron Source and Straight Guide 

The HFBS is located at the end of a dedicated straight 
neutron guide, labelled NG-2, that directly views a liq- 
uid hydrogen cold source. The cold source exhibits a 
Maxwellian spectrum with an effective temperature of 
37 K for neutrons with wavelengths 4 A < A < 15 A. |l(| 
The flux at the backscattered wavelength Ao is estimated 
to be ~ 1.2 x 10 11 n (cm 2 sSr A)" 1 . Neutrons from the 
cold source reach the HFBS by traveling along an evac- 
uated glass guide 41.1m in length with a constant 15 cm 
high by 6 cm wide cross section. The interior top and 
bottom surfaces of the guide are coated with NiCTi su- 
permirrors, which have a critical angle for reflection given 
by 9 C = Q c \ = 0.044AA -1 . The interior side surfaces of 
the guide are coated with 58 Ni-equivalent supermirrors 
for which 9 C = 0.026AA- 1 . fll| 

An 87 cm long gap interrupts the straight guide 26.3 m 
downstream from the cold source to provide space for 
filter material and a velocity selector. These elements 
are necessary because the HFBS design places the sam- 
ple position and detectors in close proximity to a direct 
line of sight with the reactor core. The estimated loss 
of Ao neutrons from this gap is « 15%. Three blocks 
of vacuum-cast polycrystalline beryllium and one block 
of "pseudo" single-crystal bismuth, each 10 cm in length 
and 16.5 cm x 6.4 cm in cross section, are used to re- 
move fast neutrons and suppress core gamma-ray (7- 
ray) radiation. [ p"2"| Vacuum-cast beryllium is known to 
produce half the beam broadening for a given length of 
filter at a given wavelength than does the hot-pressed 
grade of beryllium, filfl This implies fewer neutrons will 
be lost through the guide walls downstream after they 
have passed through the filter. The filters are also cooled 
within a liquid nitrogen dewar to minimize the effects of 
thermal diffuse scattering and maximize the transmission 
of cold neutrons. Using this 40 cm combination of beryl- 
lium and bismuth filter material, the respective transmis- 
sion of core 7-rays {E > 2MeV), and fast (E > 2MeV), 
epithermal (E ~ leV), and cold neutrons are estimated 
to be 0.20%, 0.06%, 2 x 10~ 9 %, and 53%. (The estimate 
for the cold neutron transmission includes the gap loss of 
15%.) 

The neutron velocity selector, placed just after the two 
filters, is used to limit AA/A to 18%. This minimizes the 
background produced by neutrons having wavelengths 
that lie outside the bandwidth accepted by the PST 
chopper. The blades of the velocity selector have an ax- 
ial length of 30 cm, and are composed of 0.4 mm thick 
carbon-fiber in epoxy loaded with 10 B as absorber mate- 
rial. This gives a relative suppression of unwanted neu- 
trons of 2 x 10~ 4 . The velocity selector was designed 
to accept the entire 15 cm x 6 cm neutron guide cross- 
section. This posed a technical challenge because the 
blades had to have a large enough diameter to accom- 
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modate the 15 cm vertical dimension of the guide, yet be 
able to move at high speeds. To achieve a peak neutron 
transmission centered at Ao, the velocity selector must 
rotate at 16200 rpm, which corresponds to a tangential 
speed of 410 m/s at the edge of the blades. The peak 
transmission of the selector at this wavelength is 83%. 

A standard gold foil activation analysis was performed 
along the key points of the guide with the NIST reac- 
tor operating at full power (20 MW). This analysis as- 
sumes a cross section of 98.65 barns for thermal neu- 
trons with wavelength A = 1.8 A. Just upstream of the 
filters, the thermal capture flux equivalent was measured 
to be 3.37 x 10 9 cm 2 /sec. @| After the filters and the 
velocity selector, the thermal capture flux is 1.57 x 10 9 
cm 2 /sec, which corresponds to a net transmission 
of about 44%. This value is consistent with the losses 
estimated for the guide cut and the filter and velocity 
selector transmissions. 



The performance of the converging guide is character- 
ized by the uniformity of the beam profile exiting the 
guide, and by the flux gain. To quantify these two pa- 
rameters, autoradiograph images of the beam before and 
after guide were obtained by irradiating (independently) 
two separate ~125/im thick Dy foils (1 fim = 10~ 6 m). 
The beam intensity was then integrated over 9 (5) differ- 
ent circular regions on the autoradiograph taken before 
(after) the guide, with the regions in each case represent- 
ing slightly less than 2% of the total beam cross sectional 
area. These intensities vary by 7% or less, indicating a 
highly uniform beam profile at both guide positions, j^] 
The measurements also indicate a flux gain of 3.43, some- 
what higher than the value of 3.1 predicted by the ac- 
ceptance diagram calculations. By contrast, a larger flux 
gain of 3.33 is predicted by the Monte Carlo simulations 
which take into account neutrons lost in the 87 cm guide 



B. Converging Guide 



A converging guide is used to focus the large guide 
beam cross section down to a size that is commensu- 
rate with typical sample dimensions. The guide entrance 
is located 41.3 m downstream from the cold source, just 
after the local beam shutter, and before the PST chop- 
per. All four of its interior surfaces are coated with 26^- 
equivalent (Q c = 0.044 A -1 ) supermirrors, which com- 
press the beam cross section from 15 cm x 6 cm down to 
2.8 cm x 2.8 cm. Acceptance diagrams |l5|, [l6| were used 
to determine the optimal guide length and distances over 
which the vertical and horizontal focusing should occur. 
In this context "optimal" does not necessarily imply the 
largest gain factor (although the resulting gain is still 
quite good). Instead, our design goal was to maximize 
the gain subject to the constraint that the resulting phase 
space elements corresponding to different reflections re- 
main "compact," or bunched tightly together. In so do- 
ing one avoids the presence of gaps in the divergence of 
the subsequent neutron beam seen by the PST chopper. 

Based on these considerations, the vertical focusing 
was designed to take place over a 4 m guide length, while 
the focusing in the horizontal direction occurs over only 
the last 3 m. The resulting acceptance diagrams are 
shown in Fig. || for A = Ao assuming a uniform illumina- 
tion of the guide entrance, and a perfect reflectivity for 
the supermirror coatings for < Q < Q c . Elements that 
correspond to zero and one reflection from the converg- 
ing guide walls are represented by the solid and hatched 
regions, respectively. The tiny, isolated wedge-shaped 
elements that appear in the horizontal case arise from 
neutrons that undergo two reflections. These acceptance 
diagrams were calculated analytically. As a cross check, 
Monte Carlo simulations of the straight guide + converg- 
ing guide system were performed, and these yielded ac- 
ceptance diagrams that agreed extremely well with those 
shown in Fig. |^. 
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FIG. 3: Acceptance diagrams for Ao neutrons exiting the 
converging guide. Top panel: horizontal divergence /3 versus 
lateral guide dimension. Bottom panel: vertical divergence 7 
versus vertical guide dimension. 
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cut which otherwise would have entered the converging 
guide, but would have been lost in the walls. 

Gold foil activation measurements were made at the 
entrance and exit of the guide to obtain a more accu- 
rate value of the flux gain, and the corresponding ther- 
mal neutron capture fluxes were 2.13 x 10 8 n (cm 2 sec) -1 
and 8.29 x 10 s n (cm 2 sec)" 1 , respectively, [g Q Both 
values are accurate to within 1% assuming an uncer- 
tainty of one standard deviation based on counting statis- 
tics. The resulting flux gain is then 3.89, higher than 
both the calculated value and that obtained from the 
autoradiograph images. This discrepancy motivated, in 
part, a subsequent check of the reflectivity of the straight 
guide coatings, where additional superstructure in the 
reflectivity of the 58 Ni-equivalent supermirror coatings 
was found beyond the presumed critical wave vector of 
Q c = 0.026 A -1 . Incorporating the results of the re- 
flectivity studies on the guide supermirror coatings into 
the simulations was then sufficient to reconcile the dis- 
crepancy between the experimental and calculated gains 
entirely. 

The divergence of the neutron beam exiting the con- 
verging guide has not been measured experimentally. 
However, Monte Carlo simulations of the angular dis- 
tributions of neutrons exiting the converging guide were 
done that took into account the non-zero bandwidth 
transmitted by the velocity selector. Figure ^ shows the 
results of these simulations, where N(8) is the total num- 
ber of neutrons exiting the guide with a horizontal (verti- 
cal) divergence 9 with respect to the beam axis. The ver- 
tical axes have been scaled such that the two integrated 
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FIG. 4: Monte Carlo simulations of the angular distributions 
of neutrons exiting the converging guide in the horizontal and 
vertical directions. (Simulations are courtesy of J. C. Cook.) 



distributions are equal (i. e. contain the same number 
of neutrons), but otherwise the units are arbitrary. The 
distributions at the exit of the guide are found to be 
reasonably uniform in both the horizontal and vertical 
directions. 



C. Phase Space Transformation Chopper 

Backscattering spectroscopy is inherently a flux- 
limited technique because of the narrow energy resolu- 
tion it provides. To help boost the low count rates, most 
spectrometers employ a rather poor Q-resolution that is 
introduced by the focusing analyzer system, which com- 
presses large regions of solid angle into few detectors. 
However, all cold neutron backscattering instruments at 
steady state sources are located on neutron guides, the 
coatings of which limit the divergence of the beam fed to 
the monochromator system. This then creates a situation 
that effectively defeats the purpose of the focusing ana- 
lyzer in that the monochromator system cannot supply 
a beam with sufficient divergence to take full advantage 
of the large angular acceptance of the analyzer system. 
Thus the opportunity exists to increase backscattering 
count rates, without adversely affecting the performance 
of the instrument, if one can design a primary spectrom- 
eter capable of generating a highly divergent beam from 
that supplied by the guide. To accomplish this, Schelten 
and Alefeld proposed the idea of neutron phase space 
transformation (PST) using moving mosaic crystals. flS[ ] 

The phase space transformation process is outlined 
schematically in Fig. |^. Panel (a) depicts an incoming, 
well-collimated "white" beam, similar to that which ex- 
its a neutron guide, as an element in phase space. Each 
point in this element corresponds to a neutron with an 
incident wave vector k% measured relative to the origin. 
After diffracting from a stationary mosaic crystal, the 
phase space element is transformed into the concave- 
shaped element shown on the right. No neutron energies 
are changed in this process, so there is no change in the 
number of neutrons in a given wavelength band. The sit- 
uation changes, however, if the crystal is set in motion 
along a direction perpendicular to the average scatter- 
ing vector, and antiparallel to the projection of fej onto 
the crystal Bragg planes as shown in panel (b). In this 
case the concave element rotates in phase space, and this 
does change the neutron energy distribution. More im- 
portantly, the rotation is such that shorter wave vectors 
become elongated while longer wave vectors are short- 
ened, thereby "bunching" up the wave vectors of the 
diffracted neutrons about the desired backscattered value 
ko = 1.00 A -1 . Thus the number of neutrons in a given 
wavelength band can be increased via phase space trans- 
formation. It is possible to carry this process too far by 
moving the mosaic crystals too quickly. In this case the 
outgoing phase space element rotates too much, and the 
desired "bunching" effect reverses (Fig. [5](c)). 

It can be difficult to understand how the phase space 
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FIG. 5: Simulations of the phase space transformation pro- 
cess using 0.5 cm thick HOPG crystals with an intrinsic mo- 
saic spread of 5° FWHM. The different panels correspond to 
tangential crystal speeds of (a) Om/s, (b) 250 m/s, and (c) 
500 m/s. The simulations are three-dimensional projections 
onto the (k x ,k y ) plane. 



transformation process works. An intuitive explanation 
can be given in real space keeping in mind the mosaic na- 
ture of the crystalline Bragg planes, as well as the tight 
divergence of the incident beam. The slower moving neu- 
trons incident on the PST chopper HOPG crystals must 
find mosaic blocks oriented at higher angles in order to 
satisfy the Bragg condition. Therefore they get a "push" 
from the moving crystal. By contrast, the faster neu- 
trons satisfy the Bragg condition at smaller angles, so 
diffraction occurs from crystallites oriented in the oppo- 
site sense, thereby reducing their speed. The net effect of 



this process is the remarkable conversion of the incoming 
spread in neutron wave vector (i. e. energy) into a corre- 
sponding outgoing angular spread, precisely the require- 
ment for increasing the count rate of a backscattering 
spectrometer. This transformation process is necessarily 
consistent with Liouville's theorem which requires that 
the volumes of the incoming and outgoing phase space 
elements be identical. 

Because the HFBS monochromates and energy- 
analyzes neutrons using the (111) reflection of Si, the 
PST chopper has been designed to enhance the neutron 
flux at the corresponding backscattered energy E$. As 
implemented on the HFBS, the device consists of a lm 
diameter disk whose periphery is divided into six sectors. 
Alternate sectors are covered with crystals of highly- 
oriented pyrolytic graphite (HOPG) that are 34.5 mm 
tall, 1.5 mm thick, and mounted inside protective cas- 
settes. These cassettes clamp the crystals firmly in place 
to prevent any movement which would certainly damage 
the crystals given the high speeds they experience. Fi- 
nally, the chopper disk is mounted such that the axis of 
rotation is parallel to the average graphite scattering vec- 
tor. Figure || is a photograph of the PST chopper inside 
its casing that shows the three cassettes which contain 
and protect the 180 HOPG crystals. The design of the 
PST chopper obviously introduces a pulsed-structure to 
the diffracted neutron beam, a fact that is exploited to 
good a dvan tage by the data acquisition system described 
in Sec. 



mi. 



Graphite is an ideal choice of crystal for the PST chop- 
per because the d-spacing of the (002) planes is 3.355 A, 
which is slightly larger than the 3.135 A d-spacing of the 
Si (111) reflection used to monochromate the neutrons. 
This means that neutrons that satisfy the backscatter- 
ing condition will be diffracted from the PST chopper at 
a Bragg angle of 69.2°, which allows for the convenient 
placement of the monochromator relative to the converg- 
ing guide (see Fig. ||). Moreover high-quality HOPG crys- 
tals having the requisite mosaic are readily available, and 




FIG. 6: Photograph of the PST chopper with the front cas- 
ing removed. The chopper disk and the three cassettes which 
contain the HOPG crystals are exposed. The distance be- 
tween the axis of rotation and the center of the crystals is 
505 mm. 
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the HOPG neutron reflectivity is quite good. Based 
on results from the analytical calculations and Monte 
Carlo simulations reported in Appendix |A|, the effective 
FWHM mosaic of the graphite was chosen to be 7.5°. 
This mosaic was obtained by sandwiching three crystals, 
each having a mosaic between 2.25° and 3.00°, between 
wedge-shaped spacers. This approach has the advantage 
of limiting the vertical mosaic to that of an individual 
crystal, or about 2.5°, and hence the vertical divergence. 
This is an important consideration from a design stand- 
point for two reasons. First, the vertical divergence de- 
termines the height of the monochromator. Second, given 
the high a ccelerations reached by the Doppler drive (see 
Sec. HIE), the mass of the monochromator must be kept 
a low as possible, and a larger height will imply a larger 
mass. 

A linear crystal speed of 250 m/s was chosen, again 
based on the calculations described in Appendix [A|. At 
this speed the analytical and Monte Carlo results pre- 
dict gains of 6.7 and 5.0, respectively, given the inci- 
dent neutron distributions produced by the guide sys- 
tem described earlier. To achieve a speed of 250 m/s, 
the PST chopper must rotate at 4,730 rpm. Note that 
the choice of chopper radius (0.505 m), rotational fre- 
quency (79 Hz), and number of cassettes all depend on 
the speed of Ao neutrons (630.8 m/s) and the nominal 
distance between the monochromator and the PST chop- 
per (2 m). This is because of the critical timing issue 
whereby the HOPG crystals must rotate fully out of the 
monochromated beam path by the time the neutron pulse 
diffracted from the PST chopper and then backscattered 
from the monochromator returns to the PST. With the 
aforementioned parameters, the PST chopper rotates by 
approximately 180° in this time, thereby allowing the 
now monochromatic neutron beam to pass through an 
open segment in the chopper on its way to the sample. 
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FIG. 7: Relative flux gain as a function of the tangential 
speed of the HOPG crystals in the PST. At the operating 
speed of 250 m/s (4730 rpm), the PST chopper gives a maxi- 
mal increase in neutron flux of 4.2 compared to the chopper 
at rest. 
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FIG. 8: Neutron beam intensity as a function of energy 
transfer, measured with the incident beam monitor for two 
different operating speeds of the PST. 



Because of this timing constraint, the PST chop- 
per fully transmits the monochromated neutron pulses 
at several discrete frequencies only. In addition to 
4,730 rpm, the HFBS can function with the PST chop- 
per rotating at l/3rd speed, or 1,577 rpm. To avoid these 
timing issues during the initial performance tests of the 
PST chopper, measurements of the relative neutron flux 
as a function of crystal speed were made by placing a 
detector directly on top of the PST chopper casing, i. e. 
slightly out of the backscattering condition. These mea- 
surements were performed in the fall of 1997, and are 
shown in Fig. ^. The measured gain at the operating 
velocity of 250 m/s is 4.2, somewhat reduced from the 
calculations, but still quite substantial. The HFBS is the 
first neutron spectrometer to incorporate a PST chopper. 

A very important finding of the PST chopper simu- 
lations is that the diffracted energy bandwidth has a 
FWHM of roughly 80/icV. Thus large energy transfers 
AE will come at the price of a reduced flux on the sam- 
ple, particularly at the extremes of the dynamic range. 
Because the simulations show that the distribution of 
energies reflected from the PST chopper is somewhat 
skewed to energies greater than Eq — 2.08 meV, this re- 
duction should be more pronounced for incident energies 
less than E . Fig. || displays the flux as a function of 
energy transfer AE = Ei — E f as recorded by the inci- 
dent beam monitor which is mounted between the PST 
chopper and the sample. A similar spectrum is collected 
for every data file and is used to normalize the scattered 
intensity to the incident flux. As expected from the sim- 
ulations, one sees that the number of neutrons incident 
on the sample decreases as the incident energy deviates 
from 2.08 meV. Moreover, in qualitative agreement with 
the simulation, this reduction is substantially more se- 
vere for incident energies less than E compared to those 
greater than 2.08 meV. Further details of the phase space 
transformation process are provided in Appendix |A|. 
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D. Monochromator 

Figure || shows a photograph of the HFBS monochro- 
mator which is 52 cm wide by 28 cm tall, and spheri- 
cally curved to a radius of 2.12 m. This radius of cur- 
vature is required to focus the neutrons diffracted by 
the PST chopper (2 m away) onto the sample which sits 
2.25 m from the monochromator. The monochromator 
is composed of ~ 15 hexagonal Si {111} wafers, each 
750 /im thick, glued onto the concave surface of a sup- 
port structure. Details re gardi ng the choice of crystal 
thickness are given in Sec. 1IIF which discusses the de- 
sign of the analyzer. The support structure is made of 
a graphite composite with a foam core to minimize the 
total mass (0.74 kg without crystals, 0.95 kg with crystals 
and support shaft) . It is designed to deflect no more than 
0.25 mm at the highest Doppler drive speed because large 
deflections would broaden the energy resolution. The 
two-component epoxy used to attach the silicon wafers 
to the support was selected for its durability under dy- 
namic loading conditions. The wafers were originally 
glued onto the monochromator support using a single- 
component anaerobic retaining adhesive widely used in 
the automotive industry. However, during the initial 
commissioning phase of the spectrometer there were nu- 
merous adhesive failures, particularly at high monochro- 
mator speeds. These resulted in the catastrophic loss of 
many wafers and prompted the search for an alternative. 
The two-component epoxy has proven to be very reliable, 
and there have been no failures since its first use. 

A large monochromator width is needed to span most 
of the neutron beam diffracted from the PST chopper 
which, as a consequence of the phase space transfor- 
mation process, has a much larger horizontal divergence 
(~ 15°) compared to that of the beam exiting the con- 
verging guide (~ 2.8°). The width of the monochroma- 
tor was specifically chosen so that it would capture the 
FWHM of the neutron beam diffracted from the PST 
chopper. Assuming a Gaussian angular distribution for 
the diffracted beam, this width is sufficient to intercept 
nearly 80% of the neutrons arriving at the monochroma- 
tor position. By comparison, to intercept 95% of the 
beam would require a monochromator width of more 
than 85 cm, which is impractical due to the mass limi- 
tations imposed by the Doppler drive. 

Measurements were performed to verify that most of 
the neutrons arriving from the PST chopper do indeed 
strike the monochromator. These were done by sequen- 
tially masking all but the left, middle, and right thirds 
of the monochromator surface with a neutron absorbing 
material, and then measuring the neutron intensity scat- 
tered by a standard vanadium sample with the Doppler 
drive operating at low speed. These measurements indi- 
cate that more than 40% of the monochromatic neutrons 
that scatter from the sample come from center third of 
the monochromator. This is consistent with the value of 
40% expected for a Gaussian distribution truncated at 
half-maximum. Measurements of the neutron intensity 




FIG. 9: Photograph of the HFBS monochromator, which 
is composed of Si {111} wafers glued to a light-weight, but 
rigid, graphite composite structure with a foam core. 



produced by the left and right thirds of the monochro- 
mator are symmetric, and also consistent with what is 
expected for a Gaussian distribution (~ 30% on either 
side) . 

Additional measurements were made with the 
monochromator divided vertically into thirds. These re- 
sults indicate that roughly 50% of the neutrons that 
strike the sample come from the center third of the 
monochromator, while approximately 25% come from 
each of the thirds above and below the center. If one 
assumes a vertical Gaussian distribution (which is con- 
sistent with these results) , then these measurements sug- 
gest that the monochromator is tall enough to accept 
more than 95% of the neutrons leaving the PST chop- 
per. When taken together, the vertical and horizontal 
measurements indicate that the monochromator is large 
enough to accept ~ 75% of the neutrons that leave the 
PST chopper. In fact, since the simulations described 
in Appendix A indicate that the angular distribution of 
neutrons is truncated, we believe that the fraction of 
neutrons that impinge on the monochromator is actually 
somewhat greater than 75%. 



E. Doppler Drive 

The HFBS employs a mechanical Doppler drive to pro- 
duce an oscillatory motion of the monochromator that is 
used to doppler shift the neutron energies incident on 
the sample about the average value Eg. The motion of 
the monochromator is oriented along the average silicon 
wafer [111] direction. In this manner the backscattering 
condition is always maintained, thereby preserving the 
desired sub-/ieV energy resolution. Many authors have 
discussed Bragg diffraction from moving lattices in de- 
tail. [^0|, H^] Backscattering from a moving monochroma- 
tor is a special case of this for which the neutron velocity 
is parallel to the motion of the crystal Bragg planes. In 
this case, the energy shift of backscattered neutrons rel- 



FIG. 10: Internal schematic view of monochromator doppler 
drive and cam system. From Gehring and Neumann. M 
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ative to Eq is given by 

AE = E m -E = 2E (^j + E Q {^j * , (3) 

where v m is the velocity of the monochromator. 

In contrast to other backscattering spectrometers, the 
HFBS employs a cam-based Doppler drive system to vary 
the incident energy Ei (see Fig. [lC^) . The velocity profile 
is determined uniquely by the shape of the cam, which 
is machined from tool steel and then heat treated. The 
monochromator is connected through a shaft to bearings 
that follow the contours of the cam, thereby converting 
the rotational motion of the cam into an oscillatory lin- 
ear motion. A counterweight of equal mass is driven in 
opposition by the same method to maintain the dynam- 
ical balance of the Doppler drive system. The dynamic 
range of the spectrometer is set by the frequency of the 
Doppler drive for a given cam. Cams with triangular 
and sinusoidal velocity profiles have been tested on the 
HFBS. Both have certain advantages over the other. A 
triangular velocity profile is desirable because the linear 
portions weight all energy transfers equally in time. A si- 
nusoidal waveform, on the other hand, spends more time 
at the maximum speeds, thereby weighting higher en- 
ergy transfers more heavily. This can be useful when the 
signal is weaker at large energies. The more commonly 
used crank-shaft driven Doppler systems produce a close 
approximation to a sinusoidal velocity profile. ||. 

Figure [ll] shows the triangular and sinusoidal velocity 
profiles for each cam as measured using a laser vibrom- 
eter. A small polished area in the center crystal of the 
monochromator reflects the light from the laser vibrome- 
ter. The doppler-shifted light is compared to an internal 
reference signal by the vibrometer which uses the dif- 
ference to generate an extremely accurate analog output 
that is directly proportional to the instantaneous veloc- 
ity of the monochromator. The profiles in Fig. [ll] cor- 
respond to a frequency of 4.94 Hz, or roughly a range of 
±10 /xeV in energy transfer. The velocity dependence of 
the triangular cam is linear in time over almost the entire 
cycle. In addition, the triangular cam gives a 19% larger 
dynamic range at the same frequency than does the si- 
nusoidal cam. Interestingly some "ripples" are clearly 



FIG. 11: Velocity profiles of the monochromator during one 
cycle of the Doppler drive running at 4.94 Hz. The two pro- 
files, one using the triangular cam, and one using the sinu- 
soidal cam, were measured with a laser vibrometer. The sinu- 
soidal profile has been multiplied by -1 for ease of comparison. 
Compared to the sinusoidal cam the triangular cam provides 
large regions of linear velocity dependence, and a larger dy- 
namic range at a given frequency. 



evident at higher energy transfers in both velocity pro- 
files. The ripples are real, and are attributed to the mo- 
tion of the bearings in the cam that result from the non- 
zero machining tolerances. These ripples are an intrinsic 
feature of any cam-based Doppler drive and can not be 
suppressed completely. However the vibrometer signal is 
highly localized, only reflecting the motion of the center 
of the monochromator. In fact the neutrons sample the 
entire 1500 cm 2 surface area. Since the monochromator 
support is not perfectly stiff, the ripples are not coherent 
over the surface area. Thus the effects of the ripples are 
smeared out in the reflected neutrons. As will be seen in 



Sec. |IV|, they produce no detectable effects on the energy 
spectra. 

The Doppler drive is designed to achieve a top fre- 
quency of 25 Hz corresponding to monochromator accel- 
erations in excess of 100 g. This not only poses the prob- 
lem of keeping the crystals attached to the monochro- 
mator support, but it also implies significant vibrations 
of the Doppler drive at high frequencies. Measurements 
made with an accelerometer indicate a number of pro- 
nounced resonant frequencies at which the Doppler drive 
should not be operated. The vibrations in the Doppler 



Monochromator 


Dynamic 


Energy 


Frequency 


Range 


Resolution 


4.9 Hz 


±11 ^eV 


0.80 AteV 


9.0 Hz 


±20 /xeV 


0.87 /ieV 


16.2 Hz 


±36 AteV 


1.01 AteV 



TABLE I: Instrumental elastic energy resolution as a function 
of dynamic range. 
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drive create a distribution of velocities of the monochro- 
mator around E , and therefore lead to a broadening of 
the instrumental energy resolution. For example, a mea- 
surement at 21 Hz (AE — ±46.7/^eV) gives an asymmet- 
ric instrumental energy resolution function with a width 
of about 4 /ieV (FWHM). There are several operating fre- 
quencies at which this broadening is minimized as shown 
in Table § 

F. Analyzer 

Figure [jj] shows a photograph of the HFBS analyzer. 
The analyzer consists of 8 spherical "orange peel" -shaped 
sections that stand 201cm tall by 36.4cm wide (at the 
center). Together these sections give an approximately 
continuous coverage over scattering angles 29 ranging 
from 39.3° to 124.3°. For lower angles, four Debye- 
Scherrer rings cover 7.8° < \20\ < 39°. At a coverage 
of nearly 23% of An steradians, the HFBS analyzer is 
the larger than that of any other backscattering instru- 
ment. The analyzer radius of curvature is 2.05 m, instead 
of 2 m, so that backscattered neutrons are focused onto 
the detectors rather than onto the sample. This allows 
for a smaller entrance window into the detector assembly, 
and thus a reduced background, but at a cost of a small 
increase in energy resolution and a slightly asymmetric 
shape to the instrumental energy resolution function. 

An important gain in intensity is obtained by gluing 
large diameter (~ 120 mm) Si {111} wafers onto both the 
monochromator and the analyzer support structures be- 
cause bending increases the intrinsic lattice gradient of 
silicon beyond its Darwin limit. This results in more 
flux at the expense of an increased instrumental energy 
resolution. The amount by which Ad/d changes with 
bending depends on both the radius of curvature R c and 
the crystal wafer thickness t according to the expression 




FIG. 12: Photograph of the HFBS analyzer with technician 
S. Slifer kneeling to inspect the alignment. 



ically bent Si {111}, is about 0.44. {H] 

The preceding equation predicts that a thickness of 
only 150 //m is sufficient to obtain an energy resolution 
of 0.75 fieV FWHM, which satisfies the HFBS sub-fieV 
constraint. However, the effect of a significant bend- 
ing strain on the reflectivity of perfect silicon crystals 
was not known. Therefore two sets of extensive tests 
were carried out at the Institut Laue Langevin on the 
backscattering spectrometer IN16 using small analyzer 
prototypes composed of wafers with thicknesses from 250 
to 950 fim to determine the optimal wafer thickness ex- 
perimentally. The results of these measurements indi- 
cated a much weaker dependence of the energy resolution 
on the wafer thickness t. Moreover, as shown in Fig. [l^, 
it was observed that silicon thicknesses less than 700 /im 
did not fully saturate the reflectivity. Based on these 
results a value of 750 /im was chosen for the thickness 
of the silicon wafers covering both the monochromator 
and analyzer. In addition, it was found that chemically 
etching the wafers gave a highly Gaussian shape to the 
resolution function, whereas unetched wafers resulted in 
an additional strong Lorentzian component. Hence all 
wafers used to make the HFBS analyzer and monochro- 
mator were chemically etched after cutting (but not pol- 
ished), n 

That such thick wafers are needed to saturate the 
neutron reflectivity is surprising. The primary extinc- 
tion length for the (111) reflection of silicon is only 
34.2 /tm. Popovici later pointed out that for small deflec- 
tions, a wafer simply bends and the strain remains zero 
in the middle (neutral) layer, [^ij The reflectivity will 
then saturate as long as this neutral layer is thicker than 
the primary extinction length. But when the deflection 
exceeds roughly half the wafer thickness, which happens 
very soon with thin wafers, this neutral layer becomes 
strained. The wafer stretches non-uniformly on bend- 
ing and the lattice gradient changes drastically, thereby 
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FIG. 13: Plot of neutron reflectivity versus Si {111} wafer 
thickness. Open circles represent etched wafers, while solid 
circles represent unetched wafers. The solid triangles refer to 
the two IN16 monochromators. 
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altering the primary extinction length. 

The stretching should result not just in a spread of d- 
spacings, but also in a non-uniform change in the nominal 
d-spacing, which in backscattering will matter. One then 
has to increase the thickness in order to return to an un- 
strained middle surface, and thus to a good reflectivity. 
Given R c = 2.05 m, and an effective wafer "diameter" of 
order 102mm (the wafers are hexagonal), the resulting 
deflection is about 0.63 mm. Based on this argument, 
the wafer thickness should exceed 1.2 mm in order to be 
larger than twice the deflection. Our experimental tests 
show that such thick wafers are not necessary to saturate 
the reflectivity. However they do convincingly demon- 
strate the validity of the premise that the neutral layers 
in a thin crystal wafer are drastically affected by bending. 



sample insert. These are used to detect scattering from 
the smallest three Debye-Scherrer rings, and cover mo- 
mentum transfers of 0.25 < Q < 0.47 A -1 . Another lone 
pencil-style detector with a 3 Hc fill pressure of 4 bars and 
a diameter of 25.4mm is used to detect neutrons that 
backscatter from the outer-most Debye-Scherrer ring. 
These detectors are somewhat off backscattering, result- 
ing in a slightly worsened energy resolution and an asym- 
metric instrumental resolution function. The larger ac- 
tive area and shorter collimation also lead to a higher 
background. Both the pencil detectors and the low-angle 
rectangular detectors perform well in discriminating the 
neutron signal from the 7 background. The background 
from 7 radiation and electronic noise is negligible. 



G. Detector Assembly 

Neutrons backscattered from the analyzer must pass 
through the sample a second time in order to reach the 
detectors, which are located on a 10 cm radius centered 
on the sample position (see Fig. |J). A radial collimator, 
composed of 3 cm long cadmium fins sandwiched between 
each detector, sits between the sample and detectors to 
reduce the background and suppress spurious scattering 
from the sample environment. The HFBS detector as- 
sembly consists of 12 12.27 mm diameter pencil-style de- 
tectors, each with a 3 He fill pressure of 6 bars. The detec- 
tors are located at positions corresponding to scattering 
angles of 32.3° < 29 < 121.25°, or momentum transfers 
of 0.56 < Q < 1.75 A -1 respectively. Each detector has 
an angular acceptance of 7.75°, which is about half the 
horizontal divergence of the neutron beam. The PST 
chopper housing partially shadows the analyzer at large 
scattering angles resulting in a reduced count rate in the 
detectors at Q = 1.68A -1 and Q = 1.75 A" 1 . 

A picture of the sample area taken from within the 
HFBS vacuum chamber is shown in Fig. [l4|. The detector 
assembly appears on the left-side of the cylindrical sam- 
ple insert, and is shielded by a cadmium and boronated- 
aluminum frame. The assembly and frame are electri- 
cally isolated, and mounted onto a metal plate attached 
to a stand-alone post that can be rotated away from the 
sample insert to facilitate testing and repairs. The large 
circular outline of the PST chopper is visible behind the 
sample insert, and the chopper window through which 
the main guide and monochromated beams pass is out- 
lined in black. The shiny square-shaped metal piece to 
the left of the detector assembly is the end of the beam 
stop for the main guide. At right one can see a few of 
the hexagonal silicon crystals comprising the high-angle 
portion of the analyzer. The shadowing of this portion 
of the analyzer by the PST chopper is evident. 

Three additional rectangular-shaped detectors with 
3 He fill pressures of 2.5 bar and a 38 mm x 38 mm ac- 
tive area are mounted on the PST housing, but they are 
not visible in Fig. [l4| because they are located behind the 



H. Detector Electronics 

The HFBS uses a combined preamplifier-amplifier- 
discriminator (PAD) unit to process neutron-event sig- 
nals from each 3 He detector prior to being counted by 
a scaler. Several views of the PAD unit are shown in 
Fig. [TR and the manner in which these units connect to 



the detector assembly can be seen in Fig. 14, The PAD 
design is non-trivial as it must comply with a number of 
constraints. First, the PAD units have to withstand op- 
eration under vacuum because the detectors are operated 
inside the HFBS vacuum chamber, and the PAD units are 
plugged directly onto the detectors for optimum noise re- 
duction and reliability. Second, the 1/2-inch (12.27 mm) 
diameter pencil detectors are physically arranged closely 
together to optimize coverage. This requires the PAD 
units to have a very flat form factor. Third, the scaler 
being used has ECL inputs, so the PAD units need to 
produce a digital output compatible with this standard. 

Given that normal convective cooling is not present 
in a vacuum environment, the first constraint was met 
by ensuring that the power dissipation in all compo- 
nents was low enough that even the vestigial conduc- 




FIG. 14: Picture of the HFBS sample area as viewed from in- 
side the large vacuum chamber, which houses the PST chop- 
per, analyzer, and detector assembly. 
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tive cooling through the printed circuit board, plus 
whatever radiative cooling takes place, is sufficient to 
keep all component operating temperatures within their 
limits. Two commercial (Amptek) hybrid ICs were 
chosen to perform most of the PAD functions, the 
A225 preamplifier/shaping amplifier and the A206 am- 
plifier/discriminator. Both ICs have milliwatt power dis- 
sipation and are rated for vacuum operation. The digi- 
tal output buffer is based on the Telcom TC4428A dual 
inverting/noninvcrting driver, a power CMOS part that 
also shows low power dissipation even when driving heavy 
loads. 

The required flat form factor was reached through care- 
ful component selection and package design. The final 
PAD package thickness is 11.3 mm. This was achieved 
using standard through-hole components, along with 
easy fabrication and assembly, sturdiness, and adequate 
shielding. Because the detectors were mounted so closely 
together, the locking rings on the input SHV plugs had to 
be cut off. The PAD units (and therefore the SHV plugs 
which are press-fitted into them) are held in place by a 
pin that is clamped into the detector mounting plate. 

The TC4428A digital output driver produces an inher- 
ently balanced TTL- level output. The TC4428A output 
impedance is so low that the necessary ECL-level out- 
puts can be produced by the unsophisticated but effective 
method of resistive voltage dividers to the -12V supply. It 




is expected that future applications of this driver would 
use the balanced TTL output and take advantage of the 
excellent common-mode range of commercial RS485 ICs 
at the receiving end. 



I. Operational Modes and Data Acquisition 

The HFBS employs a "real time" data acquisition 
system whose primary function is to bin the neutrons 
counted in each detector according to their incident en- 
ergy. This energy is determined by the monochroma- 
tor velocity at the time the neutron was hacks cotter ed 
from the monochromator. The data acquisition system 
also measures the angular velocity and position of the 
doppler drive cam using a digital encoder, and then cal- 
culates the monochromator crystal displacement and ve- 
locity. The velocity values are stored in a memory stack, 
and then used with the displacement values to calcu- 
late the neutron flight time to the detectors. This flight 
time information is used to set stack pointer values that 
the system will later need to retrieve the correct veloc- 
ity information when the neutrons are finally counted 
by the scalers. Further details of these calculations are 



discussed in Appendix 
which lasts ~ 20 /is, 1 



During each sampling period, 
;he scaler values are latched and 
read. Using previously calculated pointer values, the in- 
cident neutron energy (i.e. monochromator velocity) is 
retrieved from the velocity stack. The system then incre- 
ments a histogram using the appropriate incident energy 
and Q values. These tasks describe the "standard" oper- 
ational mode of the HFBS. The data acquisition system 
also supports an alternative operational mode that al- 
lows the user to perform a "fixed window scan" (FWS) 
for which the doppler drive is stopped and the instru- 
ment is simply measuring the elastic intensity as a func- 
tion of Q. This mode is generally combined with scans of 
a physical variable such as sample temperature or time. 



An example of an FWS spectrum is given in Sec. IV 



FIG. 15: Three views of the HFBS PAD unit: Top - the bare 
circuit board, Middle - a side view, Bottom - the complete 
PAD unit in its casing. 



A secondary function of the HFBS data acquisition 
system is to count only those neutrons that backscatter 
from the analyzer, and to ignore the unwanted neutrons 
that scatter from the sample directly into the detectors. 
This function is handled entirely in software, and the 
timing is controlled by the PST chopper. For each sam- 
pling period, the data acquisition system measures the 
angular position of the PST chopper using a two-pole 
resolver and then compares it to a look-up table of val- 
ues to determine if the scaler counts should be kept. To 
assist in generating this look-up table an alternative op- 
erational mode, called "teach" mode, is available that 
bins the neutrons counted in each detector according to 
the angular position of the PST chopper. As shown in 
Fig. |l^, the data obtained in teach mode clearly sepa- 
rate the huge signal generated by the unwanted neutrons 
from the much weaker signal, indicated by the shaded 
regions, that results from neutrons that have been prop- 
erly energy-analyzed. Thus "teach" mode is used to in- 
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struct the data acquisition system when to keep a neu- 
tron count, and when to discard it, based on the angular 
position of the PST chopper. Figure 16 shows one full 
period of the PST chopper which, because it rotates at 
a constant rate of 79 Hz, corresponds to 12.7msec. The 
objective of teach mode is then to adjust the positions 
of the shaded regions to achieve the highest signal-to- 
background ratio possible. 

The HFBS data acquisition system is composed of a 
number of integral hardware and software components. 
The main hardware components include a 32-bit VME- 
compatible single-board computer, a 32-bit latch, a 32- 
bit high-speed multiscaler, an optical encoder, a resolver- 
to-digital encoder, and a high-performance digital signal 
processor. The CPU provides control of the slow pro- 
cesses in the VME data acquisition. The resolver-to- 
digital encoder system provides a 12-bit digital output 
corresponding to the angular position of the PST chop- 
per. The latch is used to capture and read the encoded 
angular position of the PST chopper. The scaler accumu- 
lates counts from the neutron detectors and also counts 
the number of pulses sent from the doppler drive cam 
optical encoder, which is used to calculate the angular 
velocity and position of the cam. A high speed digital 
signal processor (DSP) is used to perform the fast VME 
readouts of the latch and scaler modules, to calculate the 
monochromator velocity and displacement, to calculate 
the pointer values, and also to increment the histogram- 
ming memory. 




FIG. 16: Neutron intensity during one full cycle of the phase 
space transformation chopper with the HFBS operating in 
"teach" mode. The regions of high intensity are caused by 
monochromated neutrons that scatter from the sample di- 
rectly into the detectors without energy analysis. Only signals 
that coincide with the time profile of the analyzed neutrons 
(shaded areas) are processed by the data acquisition system 
when operating in the "standard" mode. 



J. Background Reduction 

The HFBS design places both sample position and de- 
tector assembly in close proximity to the main guide 
beam (see Fig. ||). This was done to avoid the reflec- 
tivity loss and increased beam divergence that would re- 
sult by diffracting the beam away from the guide using 
a pre-monochromator. The extra divergence is unde- 
sirable because it would reduce the effectiveness of the 
chopper. The HFBS therefore faces the prospect of hav- 
ing a higher background than those measured at other 
backscattering spectrometers, where the neutron beams 
entering the primary spectrometers are usually already 



pre-monochromated. 
neutrons per second 



Furthermore, about 6.5 x 10 9 
exit the converging guide and 
fly towards the PST chopper. When this rate is com- 
pared to the single-detector count rate of a few counts 
per second typical for a standard experiment, the need 
for extensive background reduction is obvious. 

A number of measures have been taken to reduce the 
HFBS background as much as possible. Several of these 
have already been discussed, such as the use of liquid- 
nitrogen-cooled beryllium and bismuth filters, which re- 
duce the fast neutron background component, and the 
velocity selector, which limits the wavelength bandwidth 
seen by the PST chopper. The background is further re- 
duced by continuously flushing the flight path between 
the monochromator and the PST chopper with He gas. 
This process eliminates air scattering, which for dry air is 
7.3 % per meter, and leads to an increase in the detector 
signal. A system of masks made of absorbing boronated 
aluminum (Al loaded with ~ 5% 10 B) are used to de- 
fine the entrance and exit of the PST chopper window 
so that neutrons exiting the converging guide are either 
absorbed by the masks, or hit the rotating HOPG crys- 
tals. Neutrons incident on the HOPG crystals are then 
either diffracted, or transmitted. To absorb those neu- 
trons that are transmitted, a 2 mm thick layer of boron 
in epoxy, molded to the shape of the cassettes, is posi- 
tioned immediately behind the crystals. This boron in 
epoxy layer also functions as a beam stop when one of 
the three cassettes has rotated into the beam path. This 



is quite important because, as seen in Fig. 16, the data 
acquisition system is processing signals from the detec- 
tors precisely at this time. When none of the crystal 
sections is in the beam, the incoming neutrons are ab- 
sorbed directly after the PST chopper by a beam stop 
made from a 6 Li-bearing ceramic. An additional slit sys- 
tem positioned between the chopper and sample position 
leads to a further reduction in background. Finally, the 
stainless steel housing of the PST chopper is itself com- 
pletely covered with boronated aluminum and an outer 
layer of cadmium. 

The dominant contribution to the background, in spite 
of the features listed above, comes from the leakage of 
neutrons through the chopper exit window during the 
time a crystal cassette is blocking the beam path. We 
attribute this to multiple scattering of neutrons within 
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the chopper housing. However, this contribution drops 
by a factor of 3 when the HFBS scattering chamber is 
evacuated and the chopper entrance and exit windows 
(1mm thick aluminum) are removed, which is the stan- 
dard operating mode for the HFBS. At the same time 
we observe a 22 % increase in the scattering signal. The 
HFBS scattering chamber can be evacuated with a large 
mechanical roughing pump to a vacuum of better than 
1CP 2 mbar in less than 4 hours, at which pressure the 
high voltage to the detectors can be switched on. 



IV. PERFORMANCE 

A gold foil activation measurement at the sample po- 
sition gives a flux of monochromated neutrons on the 
sample at E of 1.4 x 10 5 n (cm 2 sec) -1 . Q The beam 
size at the sample position has been measured by irradi- 
ating a 125 /im thick Dy foil. The autoradiograph image 
is consistent with a beam size of 2.9 cm x 2.9 cm, and a 
uniform beam intensity profile. 

Figure 17 shows the spectrum of an annular vana- 
dium standard with an outer diameter of 22 mm, and 
a wall thickness of 0.88 mm. This geometry yields a 10 % 
scatterer. The spectrum was taken at room tempera- 
ture with the Doppler drive operating at a frequency of 
13.5 Hz, corresponding to a dynamic range of ±30/ieV. 
The average count rate in each detector was 106 counts 
per minute. The spectrum is shown using a logarithmic 
representation to show the Gaussian nature of the vana- 
dium spectrum more clearly. For this standard we ob- 
served a signal-to-background ratio of better than 400:1. 
Measurements on materials having a lower cross section 
for absorption than that of vanadium (e.g. samples with 
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FIG. 17: Standard vanadium spectrum measured on the 
HFBS using a Doppler drive frequency of 13.5 Hz (±30 ^teV), 
and integrating over 0.62 < Q < 1.6A -1 ). The solid line 
represents a fit to a Gaussian function plus a constant back- 
ground bg + a exp[— (x / b 2 )]. The spectrum shows an almost 
Gaussian-like energy resolution with a FWHM of 0.93 ^teV. 
Extensive efforts at background reduction have led to a signal- 
to-background ratio (a / bg) of better than 400:1. 
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FIG. 18: Dynamics of glass-forming dibutylphthalate at 
1.16 A" 1 . Structural relaxation leads to a quasielastic broad- 
ening in the scattering law S(q,u>) with increasing temper- 
ature which can be described by the Fourier transform of a 
Kohlrausch stretched exponential function (solid line). The 
spectrum at 100 K - representing the instrumental resolu- 
tion - has been scaled down by a factor of 18 to match 
S{q,LU = 0;273K). 



a high hydrogen content) give signal-to-background ra- 
tios as high as 600:1. Up to Doppler-drive frequencies 
of 16.2 Hz (|A.E| < 36/zeV for the triangular cam), all 
HFBS spectra exhibit Gaussian-like energy resolution 
lineshapes. 

One of the main applications of neutron backscatter- 
ing spectroscopy is the study of slow atomic motions in 
viscous liquids. In a first experiment the structural relax- 
ation in glass-forming dibutylphthalate was measured at 
temperatures from 100 K to 273 K using a closed cycle re- 
frigerator. p5fl Liquid dibutylphthalate was encapsulated 
inside a hollow aluminum cylinder 25 mm in diameter, 
75 mm in length, having an annular thickness of 0.08 mm. 
This geometry also yields a 10 % scatterer. The spectra 
shown in Fig. [l^ have been normalized to the monitor 
and the vanadium standard, and also corrected for a flat 
background and empty-cell scattering. The measuring 
time at 273 K was 7 hours at a count rate of 76 counts 
per minute per detector (119 counts per minute per de- 
tector at 100 K). Structural relaxation, which is respon- 
sible for viscous flow, broadens the quasielastic linewidth 
with increasing temperature. Typical for glass-forming 
systems is the stretching of correlation functions over a 
wider time scale than is expected for exponential relax- 
ations. This is an example of a system for which the 
extended dynamic range of the HFBS is especially valu- 
able as it permits the study of structural relaxation in 
viscous liquids in greater detail. 

The relaxation of the Si-0 network in sodium disili- 
catc [26 was studied at temperatures up to 1600 K using 
the HFBS high temperature furnace. Due to the chemi- 



cal reactivity of liquid silicate glasses, a platinum sample 
holder had to be used. The neutron absorption cross sec- 
tion of platinum at A = Xq is 35.9 barns. This is about 
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FIG. 19: Tunneling spectrum of 2,6-Lutidine illustrating the 
large dynamic range of the HFBS (±36/ieV), and an energy 
resolution of l.OlpeV (FWHM). 



45 times larger than that of aluminum. This large cross 
section is even more problematic given that the beam 
passes through the sample twice on its way to the de- 
tectors. Even so, the high flux provided by the HFBS 
resulted in a count rate of some 13 counts per minute 
per detector, still enough for quantitative data analysis. 

A recent measurement of the tunnel splitting in 2,6- 
Lutidine (C7H9N) served as a test of the large dynamic 
range of the HFBS as well as the excellent energy reso- 
lution. This system had been examined on a number of 
spectrometers including a backscattering instrument. j27|] 
Due to the large number of peaks in the spectrum at ener- 
gies less than 50/xeV, some of which overlap, this system 
has also been used as a test case for maximum entropy 
p8| and Bayesian analysis techniques. |^9|, The data, 
summed over 10 detectors spanning 0.62A"to 1.6 A -1 , 
and normalized to the incident beam monitor spectrum, 
are shown in Fig. fill 



As discussed in Sec. [II], the backscattering spectrom- 
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FIG. 20: A fixed- window scan of toluene that clearly shows 
supercooling. The heating and cooling rates used in this mea- 
surement were both ±0.1 K/min. 



eter can also be operated with the Doppler monochro- 
mator at rest, in which case it probes only the energy- 
resolution-limited elastic scattering intensity. The FWS 
method is a useful technique for determining under what 
conditions (e.g. temperature, time) the dynamics of the 
system being studied lie within the dynamic range of the 
spectrometer. [ ^l| In addition, this technique can be use- 
ful for probing phase transitions. As an example of this, 
the heating and cooling curves for toluene are presented 
in Fig. |2^. In this figure, the data were summed over 
10 detectors and warming/cooling rates of ±0.1 K/min 
were used. The large hysteresis on heating and cooling, 
indicative of undercooling, is clear. The solid melts at 
179 K and the liquid solidifies at about 160 K. 



V. UPGRADES AND PROSPECTS 

In March of 2002, work was completed on the instal- 
lation of a new cold source into the NCNR research re- 
actor. The complex shape of the new source, which is 
an ellipsoidal shell, allows more D2O to be introduced 
in the cryostat chamber than did the spherical shape of 
the old cold source. jl0|, |32| This modification provides 
the single biggest contribution to the overall gain in cold 
neutron flux. The center of the inner ellipsoid is also off- 
set, so that the thickness of liquid hydrogen nearest the 
source of neutrons is 30 mm, as compared to 20 mm in 
the original source. Finally, the inner ellipsoid is held 
under vacuum, whereas the inner shell of the old source 
was filled with hydrogen- vapor. The overall flux gains 
range from 40% at 2.4 A to slightly more than 100% for 
wavelengths greater than 15 A. For the HFBS, the new 
cold source has increased the usable flux at Ao by a fac- 
tor of 1.8, precisely in line with expectations. (3^] With 
this cold source upgrade, the HFBS currently provides a 
neutron flux on sample of 2.5 x 10 5 n (cm 2 sec) -1 , which 
is higher than that of any other backscattering spectrom- 
eter operating on a steady state source. 

Several ideas for improving the HFBS have been con- 
sidered and evaluated in detail. These include a new 
parabolic converging guide Q with a focal point that 
is centered on the PST chopper HOPG crystals, and a 
supermirror guide section that functions as an optical fil- 
ter. While a new parabolic converging guide is mainly 
expected to improve the signal-to-background ratio, an 
optical filter would shift the HFBS out of the direct line- 
of-sight of the neutron source, and thus allow for the re- 
moval of the beryllium and bismuth filters. The replace- 
ment of these filters with an optical filter could produce 
a substantial enhancement of the incident neutron flux 
on the sample. 



VI. SUMMARY 

The primary goal of the HFBS is to provide a sub- 
fieV energy resolution capability to the NCNR scatter- 
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ing community while maximizing the flux on sample as 
well as the dynamic range. To this end the HFBS design 
incorporates several state-of-the-art neutron optic, me- 
chanical, and electronic devices. Foremost among these 
is the phase space transformation chopper which con- 
verts the energy spread of the neutron beam exiting the 
converging guide into one that is highly divergent but 
more densely distributed about the backscattered energy 
of 2.08 meV. This device has operated flawlessly for over 
four years and enhances the neutron flux by a factor 
of 4.2. In addition, a considerable increase in signal is 
achieved using a 4 m long converging guide section and a 
large analyzer that subtends nearly 23 % of 4ir steradians. 

In contrast to other backscattering spectrometers, the 
HFBS employs a high-speed, cam-based Doppler-drive 
system to vary the incident energy Ei up to ±50 (ieV. 
The shape of the cam is machined to produce a nearly 
triangular velocity profile for the monochromator. In 
routine operation the monochromator system provides a 
Gaussian-like sub-/xeV instrumental energy resolution for 
energy transfers of up to ±36 fieV. Routine user opera- 
tion began in late 1999. Since then the HFBS has proven 
to be a reliable instrument. Up to now more than 90 
successful experiments, using temperatures ranging from 
2 K to 1600 K, have been performed by users from univer- 
sity, industrial, and research centers. At present requests 
for beam time on the HFBS exceed the available time 
by a factor of 1.8. The HFBS design and performance 
have led to sizable enhancements in both neutron flux 
and dynamic range compared to that available on other 
reactor-based backscattering instruments. 
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APPENDIX A: PHASE SPACE 
TRANSFORMATION 

The function of the PST chopper is to transform the 
shape of an incoming neutron beam in phase space such 
that the neutron flux of the outgoing diffracted beam is 
enhanced at the backscattering energy Eq — 2.08 meV. 
This idea was developed by Schelten and Alefeld to ad- 
dress the severe mismatch in divergence between primary 
and secondary spectrometers. jl8| 

Figure |^ provides a phase space diagram that outlines 
the operating principle behind the phase space transfor- 
mation chopper. The incoming neutron beam has a rela- 
tively small angular divergence that is set b y the critical 
wave vector Q c of the guide coatings, (Sec. Ill B[ ) but a 
substantial spread in wave vector Afc (k = 2ir/\). Af- 
ter diffracting from the (002) Bragg reflection of a mo- 
saic HOPG crystal at rest, the outgoing phase space el- 
ement has a much broader divergence (see Fig. ||(a)). 
Figure ||(b) demonstrates what happens when the HOPG 
crystal is set in motion parallel to — k x at a linear speed 
of 250 m/s. The outgoing phase space element is now 
both larger and more divergent. More importantly, it is 
tilted such that its arc is essentially perpendicular to kf. 
This tilt produces a much narrower energy spread about 
E . 

Analytic calculations of the gain from phase space 
transformation (PST) can be made under certain sim- 
plifying assumptions, and provide extremely useful guid- 
ance in the optimization of different physical parameters. 
The geometry for the diffraction process from the moving 
lattice is shown in Fig. |l| (reproduced from JT8|). Previ- 
ous calculations by Schelten and Alefeld were performed 
for the gain under different crystal and spectrometer con- 
figurations, but were done assuming the optimal crystal 
velocity pq| . Below we describe how the analytic equa- 
tions can be extended to arbitrary crystal velocities. 

As presented in |l8[ , the angle the diffracted phase 
space element makes with respect to the horizontal axis 
is written as 



tanx 



tan 9„ 



2 tan g tan ( 



1' 



(Al) 



where 9 g is the incident angle in the laboratory frame 
and 9 is the incident angle in the moving crystal frame. 
9 g is given by the d-spacing of the moving crystal (in this 
case graphite, hence the subscript "g") and the desired 
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, 2d, 



wavelength Ao via Bragg's law, 9 g = sin 
velocity in the crystal system is determined by 

Vk = v sin 8 g [cot 9 g — cot 9) , 



The 



(A2) 



where Vk is the velocity of the crystal and vq is the veloc- 
ity of the neutron corresponding to the desired neutron 
wavevector fcg: 



m n v = hk - 



(A3) 



Equation (A2) is obtained from the geometry in Fig. |2l] 
by noting that fcosin6* s = k* sm9 and k* cos6 = K + 
ko cos 8 g , where hK = m n VK- 



We can solve Eq. (A2) for 9 with the following result 

Vk ,_, 



tan 9 — [cot 9 g 



(A4) 



v sin t/ g 

/From Eqs. (A.4) and ( jAl| ) we can determine the tilt an- 
gle, %, of the scattered phase space element. 

/From the incident phase space element we can use the 
geometry illustrated in Fig. ^l]to relate the length of the 
element, Afc, to the quantities known so far. In terms 
of the mosaic spread, e (HWHM), (assumed known) we 
have 



tane 



Afc sin ( 



K + Ak cos 9 g + k cos 9 g 



which can be solved for Ak, 

K + kg cos 9 g 



Ak 



cot e sin 9 g — cos 9 g 



(A5) 



(A6) 



We can also find the projection of the incident phase 
space element onto the vertical axis to obtain 



Afc sin 6 g 

(K + ko cos 9 g ) sin 9 g 



sin 8 g cot e 



cos 9„ 



(A7) 



We also note that the divergence of the incident beam, 
77, and the incident beam wavevector, fc , determine the 
width of the incident phase space element 5k via 



5k = ?yfco 



(A8) 




The volume of the incident phase space element is given 
by 



V m 



2AkSk. 



(A9) 



Since the phase space volume for the incident and scat- 
tered neutrons must be equal by Liouville's theorem, we 
have 



V ps = 2Afc<5fc 
= 2Ak'5k'. 



(A10) 



/From the geometry of the scattered phase space ele- 
ment we have 



y = Ak' sin^, 



(All) 



which can be solved for Afc' since we know \ an d y from 
Eqs. (Al) and (A7) respectively, 



Afc' 



V 



smx 



(A12) 



Using Eqs. ( |A10p and (|A12[ ) we can express 5k' in terms 
of the known quantities, 



5k' 



V ps sin x 

2y 



(A13) 



A magnified view of the scattered phase space element 
is illustrated in Fig. ^| I n particular we show how the 
monochromator (placed after the PST chopper) will in- 
tersect the phase space element. Note that at the optimal 
velocity this line will be parallel to the long axis of the 
phase space element and the neutron intensity gain will 
be a maximum. From the geometry shown in Fig. |2^, 
the angle a is given by a = x + 6g + \- The projection 
of the monochromator intersection with the phase space 
element is given by min(d, 2Afc') where d is given by the 
geometry shown, 



5k' 



sin a 



(A14) 



The gain of the PST chopper is then given by the ratio of 
the final phase space width (as seen by the monochroma- 
tor) to the initial phase space width defined in Eq. 




FIG. 21: Phase space geometry for Bragg diffraction from a 
moving crystal with mosaic e (HWHM). 



FIG. 22: Magnified view of the Bragg-reflected phase space 
element. 
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FIG. 23: Relative gain in neutron flux as a function of crystal 
spee d for several different values of mosaic spread based on 
Eq. (A1S). The simulations assume a wavelength of Aq. 



and is expressed in terms of the known quantities 



gam 



min(d, 2Ak') 

•qk 



(A15) 



Fig. £3| shows the calculated gain as a function of crys- 
tal speed over a large range of mosaic spreads, 2e = 
1°, 3°, 5°, 10° and 20°. For mosaics of 3° and larger, the 
gain increases rapidly from one at low speeds, but then 
abruptly enters a linear region at a mosaic-dependent 
speed. At higher speed, the linear region abruptly ends 
and the gain decreases smoothly from this point on. For 
each mosaic value except 2e = 1°, the gain curves lie 
on top of each other at sufficiently low and high crystal 
speeds. The only difference is the location of the lin- 
ear region that clips the top of the curve. The gains 
for this linear region increase with increasing 2e, whereas 
the width of the region shrinks. The linear region occurs 
when min(d, 2Afc') = 2Ak' (see Fig. ^2|), which happens 
to be satisfied for the case 2e — 1° over the entire crystal 
speed range shown. This happens because the mosaic of 
1° is too small compared to the divergence of the incident 
beam. The maximum gain found here (19.5) occurs for 
the largest mosaic of 20° at a crystal speed of 320 m/s. 
Note however that the reflectivity of the moving crystals 
has been assumed constant for all cases of mosaic spread, 
an assumption that is incorrect for real crystals. In fact, 
at a constant thickness the reflectivity of real crystals 
decreases with increasing mosaic spread. In addition, 
these calculations assume a fc-independent, or constant 
energy spectrum incident on the PST chopper, when in 
fact the the cold source produces a Maxwell-Boltzmann 
neutron energy distribution, which is further modified 
by the guide coatings and filters downstream. These and 
other simplifications will cause the actual gain value to 
be less than that found in this analytic calculation. 

This analytical analysis quickly gives one a very good 
intuitive feel for the relative gains one can expect from 



various experimental configurations. However, as men- 
tioned above, the calculations make several approxima- 
tions, all of which tend to overstate the gain, and which 
result in the chopped-off appearance of the different gain 
curves. The most important of these approximations 
treats the scattered phase space element as a parallel- 
ogram when in fact it is a rather complicated crescent- 
shaped surface (see Fig. ||). Moreover, the locus of final 
wavevectors accepted by the monochromator collectively 
describe the surface of a sphere. Thus, instead of the gain 
being given by the intersection of a line with a parallelo- 
gram, it should be given by the intersection of these two 
more complicated curved surfaces. (The word "surface" 
is slightly misleading because both objects actually have 
non-zero width along k.) This will tend to round the 
shape of gain curves and reduce the gain, particularly at 
high speeds, due to the mismatch between the shape of 
the phase space element accepted by the monochromator 
and that diffracted by the PST chopper. To account for 
these complex geometrical effects, as well as the reflec- 
tivity dependence on mosaic spread and the true energy 
spectrum seen by the PST chopper, Monte Carlo simu- 
lations of the phase space transformation chopper have 
been performed. 

The PST chopper is located im media tely after the con- 
verging guide described in Sec. IIIB, the supermirror 



coatings of which have a critical angle of about twice that 
of 58 Ni. Based on simulations, the divergence after this 
element was taken to be 2.8° in the horizontal plane and 
4.4° in the vertical direction. The HOPG crystals were 
chosen to have a thickness of 5 mm. The incident neu- 
tron energy spectrum was chosen to be a 65 K Maxwell- 
Boltzmann distribution, in accord with measurements of 
the flux from the NIST cold source, but truncated at 
4 A to simulate the effect of a beryllium filter in the in- 
cident beam. The spectrum was also truncated at 10 A 
because wavelengths longer than this have no probabil- 
ity of being diffracted by the moving crystal. This dis- 
tribution was then multiplied by the square of the in- 
cident wavelength in order to account for the fact that 
the critical angle for total external reflection is propor- 
tional to the wavelength. The simulation included both 
horizontal and vertical mosaics, and the velocity of the 
graphite crystal. The program randomly chose the in- 
cident wavelength, the angles of incidence (within limits 
given by the divergences stated above), and the vertical 
orientation of the crystallite involved in the Bragg re- 
flection. These parameters were sufficient to determine 
the horizontal orientation required for diffraction to oc- 
cur. In addition the reflectivity of the graphite was in- 
cluded via the Bacon-Lowde equation for diffraction from 
ideally imperfect crystals. [p3| The HOPG crystals used 
in the PST chopper are most likely not all "ideally im- 
perfect." Thus even this more realistic assumption still 
overestimates the reflectivity of the HOPG crystals, and 
produces simulated gains larger than what are actually 
observed. All of these simulations have been performed 
using silicon (111) crystals as the monochromator. 
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FIG. 24: Simulations of the relative gain in neutron flux 
as a function of crystal speed for several different values of 
mosaic spread (numbers are FWHM). The simulations assume 
a 5 mm thick crystal and a wavelength of Ao- 



Two-dimensional projections of simulated Bragg reflec- 
tions from an HOPG crystal having an isotropic 10° mo- 
saic are shown in Fig. [| for three different crystal speeds. 
A coordinate system is used in which the (x, y) plane 
coincides with the horizontal scattering plane, and the y- 
axis is oriented antiparallel to the (002) scattering vector. 
The incident and final k x and k y values of the diffracted 
neutrons are represented by individual dots, while the 
value of k z is neglected. This results in dots whose total 
k vectors seem abnormally short. The reference values 
of k x and k y are indicated by the solid lines. Two effects 
are evident in these simulations. The first is that the 
phase space element volume increases as the crystal ve- 
locity (oriented opposite k x ) increases, in agreement with 
the analytic calculations. This is because the Bragg re- 
flection takes place at a lower angle in the Doppler frame. 
The second is that the diffracted element tilts in phase 
space as the crystal velocity is changed. This tilt can be 
optimized to maximize the number of neutrons that have 
the correct energy to be backscattered from a Si {111} 
crystal. 

The simulated flux gain, calculated using the same pa- 
rameters as the analytical gain shown in Fig. ^3[ is shown 
in Fig. As in the analytical calculation, the peak in- 
tensity (relative to that obtained for a crystal velocity of 
zero) is shown as a function of speed for mosaics 2e = 1°, 
3°, 5°, 10°, and 20°. There are a number of similarities 
between the simulation and the analytical calculation. 
For mosaics of 3° or larger, the relative intensity increases 
from about 1 to a broad maximum, before decreasing 
again. For the parameters chosen here, the maximum 
gain (~ 6) occurs for a crystal with a 10° mosaic moving 
at about 300 m/s. The results for a 20° crystal show a 
smaller gain due to a lower reflectivity. As in the analyt- 
ical calculation, the gain for the 1° mosaic case increases 
linearly with crystal speed, and has a relative value of 
less than one for zero speed. 



Phase-space transformation is an ideally suited method 
for boosting the count rates on backscattering instru- 
ments because it increases the divergence of the diffracted 
beam at the same time, thus alleviating the mismatch 
in angular acceptance between the monochromator and 
analyzer systems. However, there are several considera- 
tions that limit how far one can take the PST process. 
For example, a PST chopper using 10° mosaic crystals, 
which in our simulations resulted in the highest gain, 
would produce a beam with a horizontal divergence that 
is slightly larger than the 20° expected. To accept and 
use all of this divergence would require a monochroma- 
tor approximately 75 cm long. This would imply a heavy 
monochromator assembly, and present serious technical 
challenges related to the Doppler drive system. Thus 
we compromised on a more reasonable monochromator 
length of 52cm and a mosaic of about 7.5°. We also 
reduced the operational velocity of the PST chopper to 
250 m/s to make the forces on the composite disk more 
manageable and increase the reliability. Given these com- 
promises, our simulations predict that the gain from the 
PST should be about a factor of 5. This can be compared 
to the gain of 6.7 calculated in the analytical formulation 
using the same parameters, and agrees quite well with the 
experimentally measured value of 4.2. 



APPENDIX B: FLIGHT TIME OFFSET 
CORRECTIONS 

In this appendix we discuss the flight-time offset cor- 
rections (FTOSC). In particular we show that when the 
amplitude and velocity of the monochromator are appre- 
ciable, parts of the spectrum undergo an apparent shift. 
It turns out that the elastic lineshape is affected the most 
whereas scattering features at the largest energy transfers 
are affected only marginally. 

A schematic of the relevant instrument geometry is 
presented in Fig. |25|. The velocity of the neutrons 
reflected from the Doppler monochromator is a time- 
dependent quantity, v n (t) . The velocity of the monochro- 
mator, v(t), modulates the reflected velocities about the 
Bragg velocity, Vq, so that the neutrons that strike the 
sample have a velocity v n (t) = vq +v{t). For a sinusoidal 
monochromator motion with an amplitude of 4.5 cm 
and a frequency of 25 Hz, the maximum monochroma- 
tor speed is 7.07 m/s. Thus the neutron velocity varies 
between 632.7 m/s and 637.9 m/s. 

The arrival time in the detectors, i<j, of the neutrons 
that scattered from the monochromator at time t is given 
by 



t d = t + 



x(t) 



Vo Vq + v(t) 



(Bl) 



where the sample-analyzer-detector distance is d sa d — 
d sa + dad, d ms is the sample- monochromator dis- 
tance, and x(t) is the monochromator displacement (see 
Fig. |25]). If we were to neglect the motion of the 
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monochromator in calculating the arrival time (x(t) 
. v(t) = 0) we obtain the zero-monochromator motion 
arrival time, 

h=t+ dsad + drns . (B2) 

The actual arrival time, td, can be rewritten in terms of 
to as follows, 

t d = t + St, (B3) 

where 

St = dms ^ ~/° ■ (B4) 
u ° 1 + ^ 

This term, St, is the flight-time offset correction 
term. Note that it only becomes important for 
large values of the quantities x/d ms and v/vq. The 
large-amplitudc/high- frequency operation of the HFBS 
Doppler monochromator yields non-negligible values for 
St. To appreciate the magnitude of St under routine 
operating conditions on HFBS, let us assume that the 
monochromator frequency is 17.7 Hz (corresponding to a 
dynamic range of ±33/ieU). At maximum monochroma- 
tor velocity, v = 5 m/s and x = 0. At zero velocity, 
x = 0.045 m. Using d ms = 2.25 m, v = 630.8 m/s, 
and d sa d — 4.12 m, we obtain St = — 28 /is at maxi- 
mum monochromator velocity, while St = 71 /is at zero 
monochromator velocity. Since the cam that moves the 




FIG. 25: Schematic diagram of the relevant instrument ge- 
ometry which provides the definition of terms used in the 
FTOSC. 



monochromator rotates at a uniform rate, constant cam 
intervals are also constant time intervals. With constant 
time bins of 20 /is (the smallest time bin in the HFBS data 
acquisition configuration), parts of the spectrum can be 
shifted by as much as three channels. However, since St 
varies, different portions of the spectrum will be affected 
differently. 

One is usually interested in a spectrum in energy, not 
time-of-arrival. Therefore we can also calculate the im- 
plications for the energy when the FTOSC time is used 
in the data acquisition and when it is ignored. When 
the motion of the monochromator is ignored in the time- 
of-arrival in the detectors, the time-dependence of the 
energy of the neutron after reflection from the monochro- 
mator is given simply by 



E noftosc = ( B5 ) 



If we take into account the FTOSC term then 



El tosc = E(t + St). (B6) 



To first order we may expand this expression out in St, 



E^ = E(t ) + 5t—\ to (B7) 
= E + SE. (B8) 



Thus for small values of St, the spectrum will be shifted 
by an amount proportional to St. In our case with joop 
\ 7.7 Hz, we find SE — 0.3/ieV, or roughly one third of the 
instrumental resolution. When E is a maximum, dE/dt 
is zero, resulting in no shift. 

The spectral lineshape will also be distorted since the 
Jacobian in the nonlinear transformation is not unity, 
although this is a much weaker effect than the shift. As is 
the case with the shift, the distortion is most pronounced 
at the elastic peak position, where the monochromator 
velocity is zero, and smallest when the monochromator 
speed is maximum. 
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